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We propose an extension of the second-order Barker-Henderson perturbation theory for polydisperse hard- 
sphere multi-Morse mixture. To verify the accuracy of the theory, we compare its predictions for the limiting case 
of monodisperse system, with predictions of the very accurate reference hypernetted chain approximation. The 
theory is used to describe the liquid-gas phase behavior of the mixture with different type and different degree 
of polydispersity. In addition to the regular liquid-gas critical point, we observe the appearance of the second 
critical point induced by polydispersity. With polydispersity increase, the two critical points merge and finally 
disappear. The corresponding cloud and shadow curves are represented by the closed curves with 'liquid' and 
'gas' branches of the cloud curve almost coinciding for higher values of polydispersity. With a further increase 
of polydispersity, the cloud and shadow curves shrink and finally disappear. Our results are in agreement with 
the results of the previous studies carried out on the qualitative van der Waals level of description. 
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1. Introduction 

A vast majority of industrially important colloidal and polymeric materials are intrinsically polydis¬ 
perse, i.e., each particle in the system is unique in size, charge, shape, chain length, etc. This feature 
of colloidal and polymeric systems has a profound effect on their phase behavior and may cause the ap¬ 
pearance of new phases and new phase transitions. In addition, the phenomena associated with the phase 
behavior of polydisperse systems, such as fractionation, are also of technological relevance. Usually, the 
theoretical methods used to study polydisperse systems treat them as a mixture of an infinite number 
of components, each characterized by a continuous variable ^ distributed according to a certain distri¬ 
bution function /(f), e.g., for a polydisperse hard-sphere fluid, the hard-sphere size a is usually used as 
such a variable, i.e., f = cr. The theoretical study of the phase behavior of such fluids, using the methods 
of the modern liquid state theory, represents a nontrivial challenge (TJ- The main obstacle for theoretical 
description arises due to the fact that now one has to deal with an infinite number of equations for co¬ 
existing phases. One of the possibilities to overcome this obstacle is to resort to the so-caUed truncatable 
free energy (TFE) models and combine them with the possibility of their analytical description. TFE mod¬ 
els are approximate schemes, where the thermodynamic properties can be expressed by a finite number 
of generalized moments of the distribution function /(f). As a result, the formally infinite number of the 
phase equilibrium condition equations can be mapped onto a set of a finite number of nonlinear algebraic 
equations for these moments and solved using standard numerical methods. This route was recently un¬ 
dertaken in a number of studies, i.e., phase behavior of Yukawa and charged hard-sphere polydisperse 
mixtures were studied using the analytical solution of the mean spherical approximation (MSA) (ZJiU and 
high temperature approximation (HTA) (6/ More recently, HTA and dimer thermodynamic perturbation 
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theory (DIPT) for associating fluids in combination with polymer MSA were used to investigate the phase 
behavior of a polydisperse mixture of Yukawa chain molecules (7l[8). While the approaches based on the 
analytical solution of the MSA appear to be rather accurate, their application is restricted to the systems 
with a factorized version of Yukawa interaction, i.e., the matrix of the coefflcients describing the strength 
of interaction is factorized into the product of two vectors. On the other hand, HTA based descriptions, 
being less accurate, are much more flexible and can be applied to a much larger variety of the potential 
models ©. 

Attempting to improve the accuracy of the HTA based approaches, we present an extension and appli¬ 
cation of the second-order Barker-Henderson perturbation theory for the phase behavior description of 
polydisperse multi-Morse hard-sphere mixture. The paper is organized as follows: In sectionj^we intro¬ 
duce the model and in sectionj^we present a corresponding extension of the BH2 theory. Our numerical 
results for the phase behavior of the one-Morse version of the model are presented and discussed in sec¬ 
tion]^ while in section we collect our conclusions. In addition, we include an Appendix with explicit 
analytical expressions for thermodynamical properties of the multi-Morse hard-sphere polydisperse mix¬ 
ture in question. 


2. The model 

We consider the mixture with interparticle pair potential represented by the generalized multi-Morse 
hard-sphere potential 




oo, 


r^cT{^,C), 


( 2 . 1 ) 


where ^ is the polydispersity attribute, i.e., a continuous version of the species index, cr(f) is the hard- 
sphere diameter of the particle of species = [aCf) -i- a{^']]/2, z„ and cq are the the screening 

length and the interaction strength of the Morse potential, respectively. The form suggested for the multi- 
Morse potential (13 is similar to that used earlier for the multi-Yukawa potential (gj. This form is very 
flexible and can be used to model a large variety of realistic potentials by an appropriate choice of the 
coefficients Anmif,) and e.g., in reference (gJ it is used to mimic a polydisperse Lennard-Jones mixture. 
Here, Nm denotes the number of the Morse potential tails and M stands for the number of terms in the 
sum for one Morse tail. Note that original Morse potential consists of two terms, i.e., one is attractive 
and the other is repulsive. In our hard-sphere Morse potential repulsive term is substituted by the hard- 
sphere term. 

The mixture is characterized by the temperature T [or fi - (k^Ty^, where is the Boltzmann’s 
constant], the total number-density p, and by the distribution function [f f{^]d^ - 1]. 


3. Theory 

3.1. Barker-Henderson second-order perturbation theory 

To describe thermodynamic properties of polydisperse Morse hard-sphere mixture, we use here the 
Barker-Henderson second-order perturbation theory. According to this theory, Helmholtz free energy of 
the system A can be written as a sum of three terms: free energy of the reference system (Aref) and the 
two perturbation terms describing the contribution to the free energy due to Morse potential (Ai, A 2 ): 

A = Aref + Ai -I- A 2 = Ahs -I- Ai -I- A 2 . (3.1) 


Here, Aref = Ahs, where Ahs is the free energy of the hard-sphere fluid. The first-order term is: 

00 

df'p(Op(^') f drr2LrHSM(^,<':t-)g(HS)«,f';r), (3.2) 

0 
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where is the hard-sphere radial distribution function. For the second-order term, we used 

the macroscopic compressibility approximation (MCA): 


OO 

df J d^'p(ap (0 J drr 2 (^]^J[/HSM«,^':r)]"g(HS)(f,<';r), (3.3) 


where (dp/dp)^^ = is the isothermal compressibility of the hard-sphere reference fluid, which is 
obtained from the Carnahan-Starling equation and given by 


K 


HS 


( 1 -t?)^ 

1 -I- 477 -I- At]^ - 4 r]^ + 77^ ’ 


(3.4) 


where the packing fraction 77 is defined as 77 
sion for the potential ( 2 . 1 ^, we have 


: / d^piOct^l^). Substituting into 1 3.2 and 1 3.3 1 the expres- 


V 


/ p A^m Af 

d^ d^’piOpi^’) E E 

^ n=lm=l 




dz„ 


(3.S) 


PA 2 

~Y~ 


/ n Nm M 

df df'p(Op«') E E [AnmiOAnmi^’)] 

n=lm=l 


(T(^,^')G®^’(f,^': 2 z„)- 


dG^^^H(,C;2Zn) 

d(2z„) 


(3.6) 


where G*^^’(^,^':z„) is the Laplace transform of hard-sphere radial distribution function 


G‘»=^>(^,<';z„) = e^«‘^'«'>J drre-^«''g(HS)(^,^':r). (3.7) 

0 

Here, we use Percus-Yevick approximation for the hard-sphere radial distribution function, since the 
analytical expressions for its Laplace transform is known. All the rest thermodynamical quantities can be 
obtained using the expression for Helmholtz free energy (H) and standard thermodynamical relations, 
e.g., differentiating A with respect to the density, we get the expression for the chemical potential: 


mo = 


^ 1 

(PA. 

5p(0' 

i y J’ 


(3.8) 


and the expression for the pressure P of the system can be calculated invoking the following general 
relation: 


pP^pjd^piOpiO-^- (3.9) 

In the above expressions, Ahs and (^) are calculated using the corresponding Mansoori et al. expres¬ 
sions HI. 

Within the framework of the BH2 approach, the model in question belongs to the class of ‘truncat- 
able free energy models’, i.e., the models possessing thermodynamical properties (Helmholtz free energy, 
chemical potential, pressure) defined by a finite number of generalized moments. In this study, we have 
the following moments: 

777; = JdfplO 777/(0/(0, 777;(0 = Cr^ (3.10) 
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m'"’ = Jd^p(amj”^(a/(a. mf\£,)^a^q){Zn,cT), 

(p[Zn,CJ) = 

Zfi 

(3.11) 

Jd.fp(am|”Ha/(a. nif\£,)^a^q){2z„,a), 

q}[2z„,(T] 

= —^ (1 - 2z„cr - , 

(2z„)2i 

(3.12) 

= Jd^p(am‘"’(a/(a, m|"’(a = o'[(pi(z„,cr), 

(pi{z„,(j] = 

(3.13) 

= Jd^p(am‘"Ha/(a, m‘”’(a = o-[^i(2z„,(7), 

(pi{2z„,a) = 

(3.14) 

d^p(amp(a/(a, 

m 

(nm)(a _ 

(3.13) 

d^^p(amp(a/(a, 

fh 

(«m)(a _ 

(3.16) 


Closed form analytical expressions for thermodynamical properties (Helmholtz free energy, chemical 
potential, pressure) in terms of the generalized moments (3.10|l-l|3d^ are presented in the Appendix. 


3.2. Phase equilibrium conditions 

The main obstacle in theoretical studies of the phase behavior in polydisperse systems arises due to 
the fact that one has to deal with an infinite number of equations for coexisting phases. However, for the 
present ‘truncatable free energy model’, these equations can be written as a finite number of equations 
for the corresponding generalized moments of the distribution function /(^ flol . 

We assume that at a certain temperature T, the system, which is characterized by the parent density 
and parent-phase distribution function separates into q daughter phases with the densities 

p**?*, and q daughter distributions • • •> All phase equilibrium condi¬ 

tions of the polydisperse system can be simply obtained by generalizing from the multicomponent case 
via the prescription p,- p/(f)d^. Due to this substitution, thermodynamic properties (Helmholtz free 

energy, chemical potential, pressure) become functionals of the distribution function f(^]. Thermody¬ 
namic conditions of phase equilibrium imply the equality of the pressures, 

p™ (r, [/™ (01) = (r. (01) = ■ ■ ■ = «)]). (3.17) 

and of the chemical potentials for each species 

(^, T, [/‘i> (ai) = pP (a T, [/'2) (ai) =... = (a a (ai). ( 3 . 18 ) 

The phase separation constrained by the conservation of the total number of particles of each species a 

/('»(a = /'“’(ax^“\ (3.19) 

a=l 

where is the ratio of the total number of particles, In phase a to the total number 

of particles In the parent phase . The conservation of the total volume occupied by the parent phase: 

(3.20) 

a=l 

where = l/p^“’, (a = 1,2,3,...,< 7 ). Finally, the normalization of the /^“’(a 

J/“^(ad^=l (3.21) 

in equation (3.19| implies the conservation of the total number of particles: 

(3.22) 

a=l 
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In this study, we will consider a two-phase fractionation of polydisperse Morse hard-sphere ntixture 
(a = 1,2). We assume that thermodynamical properties of the model depend on K generalized moments 
mo, nil, ni 2 ,..., which are defined as follows: 

mfc = p J mfc(0/(Od^. k^O, (3.23) 

and mo = p. For the case of two-phase equilibrium (a = 1,2), the conditions |3.17| - (3?^ lead to the 
following set of relations: 




(3.24) 



for a = 1 or a = 2, 


(3.2S) 


and 


md'(^)/f'’^(f)H(^, r, m®,{md^}{m*°'})d^, k^O, 


(3.26) 


where 


H(^,r,m^2’,{m™}{m™}) = 


(pd^ - p^2^] Ai2{£,, T, mQ^',{md)}{m*°H) 

(p(i)p(2)/p(0) _ p(2)j + (p(i) _ p(i)p(2)/p(0)) T, m®, {mtFllmto)}) 


(3.27) 


Ai2(^,r,{md)}{m'2J}) 



exp r, {m'^H) - p™ (^, T, {m^)}) 


(3.28) 


and Pex is the value of the chemical potential in phase 1 in excess of its ideal gas value. 

The solution of the set of equations j3.24) - jT^ for a given temperature T, for the density of the par¬ 
ent phase and for the parent species distribution function gives us the coexisting densities 

p^“^ of the two daughter phases and the corresponding species distribution functions (^), a = 1,2. The 
coexisting densities for different temperatures give us binodals, which are terminated at a temperature 
for which the density of one of the phases is equal to the density p^”^ of the parent phase. These termina¬ 
tion points form the cloud and shadow coexisting curves that intersect at the critical point characterized 
by the critical temperature Ter and the critical density per = p^' = p*^^ = p*°\ The cloud and shadow 
curves can be obtained as a special solution of the general coexisting problem, when the properties of 
one phase are equal to the properties of the parent phase: assuming that the phase 2 is the cloud phase, 
i.e., p*^' = p*°\ and following the above scheme we wdl end up with the same set of equations, but with 
p^^^ and (f) substituted by p^**^ and (^), respectively. 


4. Results and discussion 


In this section we present numerical results for the phase behavior of polydisperse Morse hard-sphere 
mixture. For the species distribution function /(f), we have chosen a log-normal distribution, i.e.. 


^(LN)(^)^ 


I 


\/27ilnI 


exp 


ln2(j3/2^) 

21n7 


(4.1) 


where / is the polydispersity index. Log-normal distribution frequently occurs at the colloidal and poly¬ 
meric processing (T). Note that in the monodisperse limit (/ = 1), this distribution is represented by the 
Dirac delta-function, d(f -1). On the contrary, when 1 becomes very large (7 » 1), the above distribution 
becomes very wide, increasing hereby the importance of the particles with a large value of f. All calcu¬ 
lations were carried out for the one-Morse version of the pair potential & i.e., Nm = 1, M = 1, with 
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zi = I.Sctq. We consider two types of polydispersity of the model: polydispersity only in the strength (am¬ 
plitude) of the pair potential Ylii(^) and polydispersity in both, amplitude and hard-sphere size 

In the former case, we have chosen (0 = and cr (0 = cto. while in the latter case (f) = 
and ct( 0 = Here, zIq = 1 and cto is the hard-sphere size for a monodisperse version of the model 

at 7 = 1, which is used as a distance unit. In what follows the density p and temperature T are presented 
in reduced units, i.e., p* - pcjg and T* - k^TI cq. 

As a first step in our numerical study, we perform the calculation of thermodynamical properties 
of the monodisperse version of the model, i.e., for 7 = 1 using the present version of the second-order 
BH theory and the reference hypernetted chain approximation with the bridge function due to Verlet 
(TTl[T2l. The latter theory is known for being very accurate in predicting the properties of simple fluids 
ca. The comparison of the results of both theories for the pressure and chemical potential (figure at 
three different temperatures demonstrate that BH theory is capable of giving relatively accurate results 
at lower and intermediate densities. At higher values of the densities, the predictions of the BH approach 
are a bit less accurate. 


6 
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Q- 
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1 

0 


-1 

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
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Figure 1. The pressure as a function of the density (the upper panel) and the chemical potential as a 
function of the density (the lower panel) along three isotherms; the upper set of curves corresponds to 
T* - 2.5, the intermediate set refers to T* = 2 and the upper set belongs to T* = 1.5. Crosses are RHNC 
results, solid lines correspond to the Barker-Henderson second-order results. 

Next, we perform the calculation of the phase diagram of a polydisperse version of the model at 
different values of polydispersity index 7. Our results for the model with polydispersity only in the am¬ 
plitude of the potential are presented in figurej^and for the model with polydispersity in both amplitude 
and hard-sphere size are shown in figurej^ One can see that polydispersity has a profound effect on the 
phase diagram. When 7=1, the phase diagram consists of the usual binodal ending in the critical point. 
In this case, the cloud and shadow curves coincide, the critical point being located at their maximum. 
For a polydisperse system (7 A 1), for each parent phase density p^°\ there is a different binodal. Each 
binodal is truncated at a maximum temperature, with the corresponding densities, and p^^^ lying 
on the cloud and shadow curves, respectively. For a critical value of p*°\ p^**^ = pcr, the corresponding 
binodal passes through the intersection of the cloud and shadow curves. This occurs at T = Ter and, since 
Per - p'-^^ = p^^^ = p*°\ the point (per, 7’cr) is a critical point, where two coexisting phases, ( 1 ) and ( 2 ), 
become identical. We are interested in the phase behavior of the system at relatively large values of poly- 
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Figure 2. Phase diagrams of the polydisperse Morse hard-sphere mixture with amplitude polydispersity 
only in the ip* ,T*) plane for three different values of the polydispersity index I, I = 1.02967 (the upper 
panel), I = 1.034 (the intermediate panel) and / = 1.035 (the lower panel), obtained using the Barker- 
Henderson second-order perturbation theory, which includes cloud (soUd line) and shadow (dotted line) 
curves, two critical points and critical blnodals (dashed lines). Filled circle denotes the position of the 
critical points. 


dispersity. For a small polydispersity, the system has only one critical point |2}{6], which originates from 
the regular liquid-gas (LG) critical point of the corresponding monodisperse version of the system. With 
the polydispersity increase, there appeares an additional critical point induced by polydispersity. This 
effect on the qualitative van der Waals level of description has been observed by Bellier-Castella et al. 
(^|l3. The second critical point, which we denote as polydisperse (P) critical point, is located at larger 
values of the density and at lower values of the temperature. It is present for both tjqtes of polydisper¬ 
sity studied (figure and [^. With polydispersity increase, both LG and P critical points move towards 
each other and, for a certain limiting value of polydispersity, they merge. There are no critical points 
above this limiting value (figure|^ lower panel). With a further increase of polydispersity and at lower 
temperatures, we expect that the two-phase coexistence becomes unstable and there appeares a region 
of three-phase coexistence. For relatively high values of polydispersity studied here we observe a rather 
unusual shape for the cloud and shadow curves. For both types of polydispersity, they are represented 
by closed curves of elipsoidal- and A-like shapes for the shadow curves seen in figure and figure 
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p* 



P=K 



p* 

Figure 3. Phase diagrams of the polydisperse Morse hard-sphere mixture with size and amplitude polydis- 
persity in the (p*, T*) plane for three different values of the polydispersity index I, I = 1.035 (the upper 
panel), I = 1.045 (the intermediate panel) and / = 1.055 (the lower panel), obtained using the Barker- 
Henderson second-order perturbation theory, which includes cloud (sohd line) and shadow (dotted Une) 
curves, two critical points and critical binodals (dashed lines). Filled circle denotes the position of the 
critical points. 


respectively, and closed curves of the linear shape for the cloud curves (figure [^[^. In the latter case, 
the ‘liquid’ and ‘gas’ branches of the cloud curves almost coincide for the larger polydispersity (figure]^ 
and[^ lower panels). With a further increase of polydispersity, the cloud and shadow curves shrink and 
finally disappear. 

Finally, In figure]^ we display distrlbufion functions of the polydisperse Morse hard-sphere mixture 
with size and amplitude polydispersity at two values of the temperature, which are higher and lower than 
the second critical point temperature, i.e., T* = 1.9 and T* - 1.4, respectively. We present distribution 
functions of the coexisting phases on the critical binodal at T* -lA (figure]^ upper panel) and on the 
coexisting cloud and shadow phases at T* - 1.4 (figure intermediate panel), and T* - 1.9 (figure 
lower panel). As usual (HO, on the binodal the particles with larger values of ^ fractionate to the liquid 
phase while particles with smaller values of ^ fractionate into the gas phase. Fractionation of the particles 
to the shadow phases for the model at hand depends on the temperature. For the temperatures larger 
than the temperature of polydispersity induced critical point, we observe fractionation of the usual type, 
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Figure 4. In the upper panel, we display parent (f) (solid line) and daughter (f) (corresponds to 
gas phase, dashed line) and (f) (corresponds to liquid phase, dotted line) distribution functions for 
critical binodal at / = 1.055 and T* - 1.4. In the intermediate panel, we show (f) (corresponds to gas 
shadow phase, dashed hne), (f) = (f) (corresponds to liquid and gas cloud phases, solid line) and 

(^) (corresponds to liquid shadow phase, dotted line) distribution functions at / = 1.055 and T* = 1.4. 
In the lower panel, we show the same as in intermediate panel at T* = 1.9. 


i.e., the liquid shadow phase contains particles with ^ larger than those of the gas cloud phase while the 
gas shadow phase contains particles with ^ smaller than those of the liquid cloud phase (figure]^ lower 
panel). Note that distribution functions for the gas and liquid cloud phases are always the same and equal 
to the distribution function of the parent phase. The situation changes for the temperatures lower than 
the temperature of the second critical point (figure]^ intermediate panel). In this case, both liquid and 
gas shadow phases contain particles with lower values of ^ than those of the liquid and gas cloud phases. 
At the same time, the liquid shadow phase has particles of a larger value of £, than the gas shadow phase. 
This behavior of the model is related to the fact that for temperatures lower than the temperature of the 
second critical point, both branches of the shadow curve are located to the left of both branches of the 
cloud curve, i.e., the density of the shadow phases is always lower than the density of the cloud phases. 
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5. Conclusions 

In this paper we present an extension of the second-order Barker-Henderson (BH2) perturbation the¬ 
ory for a polydisperse hard-sphere multi-Morse mixture. To verify the accuracy of the theory, we compare 
its predictions for the limiting case of a monodisperse system, with predictions of the very accurate ref¬ 
erence hypernetted chain approximation. The theory is used to describe the liquid-gas phase behavior 
of the mixture with different type and different degree of polydispersity. In agreement with the previous 
study |Io]|T 3|, which was carried out using the qualitative van der Waals level of description, we observe 
the appearance of the second critical point induced by polydispersity. With the polydispersity increase, 
the two critical points merge and finally disappear, i.e., for polydispersity larger than a certain thresh¬ 
old value, there are no critical points. The corresponding cloud and shadow curves are represented by 
the closed curves with ‘liquid’ and ‘gas’ branches of the cloud curves that almost coincide with the poly¬ 
dispersity increase. With a further increase of polydispersity, the cloud and shadow curves shrink and 
finally disappear. 


Appendix 


Here, we present expressions for thermodynamical properties in terms of the moments (3.10^ - j3l^ . 
We have: 




zlD^^hzn) dZn 


zlD[^\zn) 


■ \ QT'' + QT'' 

zlD^^hzn) I ^ 0 


2 

-h 


dDf(Zn) 


D''^\zn) dZn 


(A.1) 




zlD[^\2Zn) d{2Zn) 


Q\^"'\2Zn) 

[2z„)^Dl"\2Zn) 




~(nm) dQ^^""‘>{2Zn) 

, Q ---t Q\"’^\2zn) 

4zlD^f{2Zn)^ ^ “ 


J_ 1 

■Zn D‘"’(2 z„) d[2Zn) 


where 


Df {z„] = a2 - (a + + ^) - 2;r I Am'"> + ^ [m 2 + +2m'"’) - 2 


, (A.2) 


(A.3) 


A= 1- 


7tm^ 




Anm).2 


[m3 -t z„m^”’j -t [y m2 -t Znm!-"^'! [m'""®^j 

+ z„(A-7 rm‘">)m‘”"*>m'"'">|, 


(A.4) 


Q^^hzn) = |Am‘”'")m'”™) + | [(m3 + z„mf) m‘”'”>m'"'”> + (|^m2 + z„m'”)) m‘”'”>m, 
,(A-7rmM)(m'”'”’)H, 


{nm) 

2 


+ Zn 


(A.S) 


dQ'^a"'\zn) n 
dZn 2 


dm. 


dZn 


(fi) 

" +m“ 


(<«)’ 


dm, 


(fi) 


- 1 Zn—- 

2 dZn 


+ m' 




dm'"’ 

^ dz, 


(A.6) 


13605-10 













































Second-order Barker-Henderson perturbation theory 
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SacTOcyeaHHJi TepMOAi/iHaiviiMHoY Teopii 36ypeHb APyroro 
nopjiAKy EapKepa-XeHAepcoHa A-nn AOC-niA>K6HH5i (|)a30B0T 
noeeAiHKi/i noniAi/icnepcHoY cyMimi reepAi/ix C(|)ep Mop3e 

T.B. rBOBflb, fO.B. Ka.nKiXHHii 

iHCTMtyT c|)i 3 i/iKM KOHfleHCOBaHi/ix ci/iCTew HAH VKpaiHn, By.fl. 1. CBeHpipbKoro, 1, 79011 /IbBiB, yKpai'Ha 

BanponoHOBaHO sacrocyBaHnn TepMOflMHaMiHHoT reopiT sBypeub flpyroro nopnflKy BapKepa-Xenflepcona fl/in 
flOcni/pKeHHn noniflncnepcHOi cyMiini TBepfli/ix cejjep Mopse. Ann nepeBipKi/i TonnocTi nopiBHioiOTbcn pe3y.fib- 
ratM piei Teopiifl/in rpaHi/iHHoro BMnaflKy MOHOflncnepcHOi ci/icreMn 3 pe3y/ibTaTaMM flyxe Tonnoro 6a3ncHoro 
rinepnaHpioxKOBOro Ha6/iM>KeHHn. Teopin Bi/iKopncTOBycTbcn fl/in oni/icy 4)a30B0i noBefliHKi/i piflHHa-ra3 fl/in cy¬ 
Mimi 3 pi3HMMM TMnaMM Ta pi3HMMM CTyneHnMM nOniflUCnepCHOCTi. OxpiM 3BMHariHOi Kpi/ITl/mHOiTOnKM pifll/IHa- 
ra3, MM cnocrepiracMO noney flpyroi kpmtmhhoi toukm, axa e ayMOBnena no.fiiflMcnepcHiCTio. I3 36mbmeHHnM 
n0.fliflMCnepCHOCTi pi flBi XpMTMHHi TOHKM 3.f1MBaiOTbCn i, Hapemti, 3HMKaiOTb. BiflnOBiflHi KpMBi XMapM Ta Tini 
npeflCTaBneni aaMKHOHMMM kpmbmmm 3 rmKaMM pifli/ina Ta raa, flnn xpMBOi XMapn bohm Marixte 36iraiOTbcn fl/in 
BMupMX 3HaHeHb noniflMcnepcHOCTi. flpM noflanbmoMy 36inbmeHHi noniflMcnepcHOCTi xpMBi XMapM Ta Tini cko- 
ponyiOTbcn i, napemTi, 3HMKaiOTb. Hami peaynbTaTM yaroflxcyiOTbcn 3 pe3y./ibTaTaMM nonepeflnix AOC/iiflxenb, 
nxi 6y.fiM npoBefleni na nxiCHOMy pieni onMcy Ban ^ep Baanbca. 

K/iioMOBi cnoBa: TepMOfli/tHaMiHHa reopifi aSypeHb, nonifli/icnepcHicTb, (paaose cnisicHyBaHHn, Kono'iflHi 
cucreMi/i, noTenpia/t Mopae 
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